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ABSTRACT 

We present Rossi X-Ray Timing Explorer and Swift observations made during the final three weeks of the 
2006-2007 outburst of the super-Eddington neutron star (NS) transient XTE J170 1-462, as well as Chandra 
and XMM-Newton observations covering the first ~800 days of the subsequent quiescent phase. The source 
transitioned quickly from active accretion to quiescence, with the luminosity dropping by over 3 orders 
of magnitude in ~13 days. The spectra obtained during quiescence exhibit both a thermal component, 
presumed to originate in emission from the NS surface, and a non-thermal component of uncertain origin, 
which has shown large and irregular variability. We interpret the observed decay of the inferred effective 
surface temperature of the NS in quiescence as the cooling of the NS crust after having been heated and 
brought out of thermal equilibrium with the core during the outburst. The interpretation of the data is 
complicated by an apparent temporary increase in temperature ~220 days into quiescence, possibly due to 
an additional spurt of accretion. We derive an exponential decay timescale of ~ 120^20 days for the inferred 
temperature (excluding observations affected by the temporary increase). This short timescale indicates a 
highly conductive NS crust. Further observations are needed to confirm whether the crust is still slowly 
cooling or has already reached thermal equilibrium with the core at a surface temperature of ~125 eV. The 
latter would imply a high equilibrium bolometric thermal luminosity of ~5 x 10 33 erg s _1 for an assumed 
distance of 8.8 kpc. 

Subject headings: accretion, accretion disks - stars: neutron - X-rays: binaries - X-rays: individual (XTE 
J1701-462) 



1. INTRODUCTION 

Understanding the internal properties of neutron stars 
(NSs) remains one of the major unresolved problems in as- 
trophysics. Efforts to constrain these properties most com- 
monly focus on narrowing down the allowed regions in NS 
mass-radius diagrams, to thereby rule out some of the pro- 
posed equations of state for the matter inside the stars. An 
alternative a pproach is to observe the cooling of NSs (for 
a review see lYakovlev k. Pethickl 120041 ) . Initially, this ap- 
proach focused mainly on the long-term cooling of isolated 
NSs over the first ~10 6 years after their birth, but it has in 
the past decade been extended to NSs reheated by transient 
accretion. 

As matter from a binary companion is accreted onto 
the surface of a NS, matter already present is compressed 
further down into the crust to higher densities. This 
leads to heating due to nuclear reactions, so-called deep 
crustal heating Brown et all 119981 iRutledge etHI l2002t 
lHaensel fc ZdunikJl2008f ). Most of the heat is produced by 
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pycnonuclear reactions hundreds of meters below the sur- 
face and is spread throughout the star by heat conduction. 
Cooling takes place via neutrino emission from the interior 
and photon emission from the surface. In ~10 4 years the 
NS enters a limit cycle in which heating during accretion 
episodes is on av erage balanced by cooling during and be- 
tween outbursts (|Colpi et alj 12001). The temperature of 
the NS core is not expected to change appreciably after 
this, and its value depends on the long-term time-averaged 
mass accretion rate as well as on the efficiency of the cooling 
mechanisms at work, which in turn depend sensitively on 
the properties of the material inside the star. High-mass 
NSs are thought to potentially have much more powerful 
neutrino emission mechanisms active in their cores, com- 
pared to their low-mass co unterparts (|Yakovlev et al.| [2003: 
lYakovlev fc Pe thick 2004); this is referred to as enhanced 
cooling, in contrast to the so-called standard cooling of the 
low-mass NSs. 

The blackbody-like component often seen in spectra from 
NS low-mass X-ray binaries (NS-LMXBs) in quiescence is 
usually interpreted as thermal radiation from the surface 
of the NS. For a NS in thermal equilibrium, the tempera- 
ture of the core can be inferred (via a model of the crust 
temperature profile) f rom the effective surface temperature 
(|Yakovlev et al.lr2004T ). A measurement of the thermal sur- 
face emission, in conjunction with information on the dis- 
tance to the source and the average length, intensity, and 
recurrence time of accretion episodes, can thus possibly be 
used to constrain the properties of the material in the core 
of the NS. 

In contrast to the core, the NS crust can have its temper- 
ature significantly altered, and be brought out of thermal 
equilibrium with the core, during single outbursts. In most 
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NS X-ray transients, outbursts last from weeks to months, 
with periods of quiescence in between lasting from months 
to decades. The temperature of the crust is then expected 
to be raised only slightly during outbursts, and thermal 
equilibrium with the core will be quickly re-established 
(within days to weeks). However, in the so-called quasi- 
persistent transients, outbursts can last for several years or 
decades, in which case the crust is expected to cool down 
to equilibrium on a m uch longer timescale of several years 
(|Rutledge et al.l 120021 ) . It is therefore possible to monitor 
the cooling of such quasi-persistent transients with satel- 
lites such as Chandra or XMM-Newton. The timescale of 
the cooling is dependent on the properties of the material in 
the crust, such as its thermal conductivity, and structures in 
the cooling curve can give information about the nature and 
locat ion of heating sources in the crust (jBrown fe Cummind 
12001 . 

Since the advent of Chandra and XMM-Newton, only 
a handful of NS transients have entered quiescence after 
long-duration (year or longer) outbursts. KS 1731-260 
and MXB 1659-29 entered quiescence in 2001 after out- 
bursts lasting around 12.5 and 2.5 years, respectively. Both 
sources were observed to co ol down to a constan t level 
over a period of a few years (f Wiinands ~et~aT1 120011 [2002; 
WiinandsH2002HRutledge et al.ll2002t ;lWiinands et al.1 12003. 



20041: IWiinandsll2004t ICadrett" et al. 20061 120081) ." though a 
recent observation of KS 1731-260 at more than 3000 days 
post-outburst suggests it may still be cooling slowly (E. M. 
Cackett et al. 2010, in preparation). The observed cool- 
ing timescales were interpreted to imply a high thermal 
conductivity for the crust, in agreement with more recent 
findings fr om the fitting of theoretical models to the coo l- 
ing curves (|Shternin et al.ll2lM iBrown fc Cummingil2009D . 
In 2008, EXO 0748-676 entered quiescence after active ac- 
cretion for over 24 years. Swift and Chandra observations 
of the source in the first half a year since the end of the 
outbu rst indicate very slow initial cooling ([Degenaar et al.l 
2009). In contrast to KS 1731-260 and MXB 1659-29, EXO 
0748-676 has shown a significant non-thermal component 
in its spectra in addition to the thermal component. Such a 
non-thermal component has been seen for many quiescent 
NS-LMXBs. It is usually well fitted with a simple power law 
of photon index 1-2 and typically domi nates the spectrum 
above a few keV (|Campana et al.lll99 8a). A number of qui- 
escent NS sources have spectra which are completely dom- 
inated by the power-law component and do not require a 
thermal comp onent, e.g., the milli second X-ray pulsar SAX 
J1808.4-3658 dHeinke et al.ll2007D and t he globular cluster 
source EXO 1745-248 ([Wiinands et aLl feOOS). The power- 
law compo nent is common amon g millisecond X-ray pulsars 
(see, e.g., iCampana et al.l [20051 ). but its origin is poorly 
understood. Suggested explanations include residual accre- 
tion, either onto the NS surface or onto the magnetosphere, 
and a shock from a pulsar wind (see, e.g.. ICampana et ah! 
Il998al) . We note that it has also been argued that low- 
level spherical accretion onto a NS surface can prod uce a 
spectrum with a thermal shape (|Zampieri et al.lll995l ). 

1.1. XTE J1701-462 

XTE J1701-462 (hereafte r J 1701) was disco vered with 
the All-Sky Monitor (ASM; iLevine et ahl 119961 ) on board 
the Ross i X-Ray Timing Explo rer (RXTE) on 2006 Jan- 
uary 18 ([Remillard fc Lin1l2006l ). shortly after entering an 
outburst (see Figure Q}. Re-analysis of earlier ASM data 
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Figure 1. RXTE ASM light curve of XTE J1701-462 showing the 
2006—2007 outburst and the subsequent quiescent period. Data points 
represent 1 day averages. The upper row of vertical bars indicates the 
times of the ten Chandra observations made after the end of the out- 
burst; the lower row indicates the times of the three XMM-Newton ob- 
servations. No other observations of XTE J1701-462 sensitive enough 
to detect the source have been made since the outburst ended. 

further constrained the start of the outburst to a date be- 
tween 2005 December 27 and 2006 January 4 (Homan et al. 
120071 ). During the -1.6 -year-long outburst the source be- 
came one of the most luminous NS-LMXBs ever seen in the 
Galaxy, reaching a peak luminosity of ~ 1.5 £Edd, an d it ac- 
creted at nea r-Eddington lumi nosities throughout most of 
the outburst (|Lin et al.ll2009bD . The source entered quies- 
cence in early 2007 August (see Section 12.61 for a discus- 
sion of our definition of quiescence for this source). Dur- 
ing the outburst the source was monitored on an almost 
daily basis with RXTE. Spectral and timing analysis of the 
early phase of the outburst is presented in iHoman et al.l 
(I2007I) . and lLin et all ((2009b) give a detailed spectral anal- 
ysis of the entire period of active accretion. In the early 
and most luminous phase of its outburst, J 1701 exhib- 
ited all spectral and timing characteristics typical of a Z 
source, and is the only transient NS-LMXB ever observed 
to do so. During the outburst the behavior of the source 
evolved throu gh all spectral subclasses of low -magnetic-field 
NS-LMXBs (jHasinger fc van der Klisl I1989T ) . starting as a 
Cyg-like Z source, then smooth ly evolving into a Sco-like 
Z source (Kuulkc rs et al.l[l997l) . and finally into an atoll 
source (first a bright GX-like one and subsequently a weaker 
bursting one). This evolution will be discussed in detail in 
an upcoming paper (J. Homan et al. 2010, in preparation). 
The unique behavior of the source in conjunction with the 
dense coverage by RXTE has made it possible to address 
long-standing questions regarding the role of mass accre- 
tion rate in causing t hese subclasses an d the spectral states 
within each subclass ([Lin et al.H2009bD . Toward the end of 
the outburst J1701 exhibited three type I X-ray bursts, the 
latter two of whic h showed c l ear pho tospheric radius expan- 
sion. From these iLin et ail (|2009a| ) derive a best-estimate 
distance to the source of 8.8 ± 1.3 kpc, using an empiri- 
cally de termined Eddington l uminosity for radius expansion 
bursts (|Kuulkers et al.ll2003f ). 

J1701 provides a special test case for NS cooling. It ac- 
creted for a shorter time than the three cooling transients 
with long-duration outbursts mentioned above, but for a 
longer time than regular transients. Moreover, the level 
at which it accreted is higher than for any other NS tran- 
sient observed. This source therefore allows new parameter 
space in NS cooling to be probed. The close monitoring of 
the source with RXTE also makes it possible to get a good 
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estimate for the total fluence of the outburst. This gives 
information about the total mass accreted and hence about 
the heat generated from crustal heating, a crucial input pa- 
rameter for theoretical models of the cooling. Flux values 
derived from spectral fits to RXTE data (spectra from 32 
s time bins, with linear interp olation between data points; 
see Fig. 3 in lLin et al.ll2009"bfi imply a total bolometric en- 
ergy output (corrected for absorption) during the outburst 
of ~1.0 x 10 46 erg for an assumed distance of 8.8 kpc and 
system i nclination of 70° (D. Lin 2009, private communica- 
tion; see lLin et al.ll2009bl for details on the spectral fitting). 
This value is likely to be uncertain by a factor of ~2-4 due 
to uncertainties in the distance and inclination of the sys- 
tem, as well as in the choice of a correct spectral model and 
the extrapolation of the model outside the RXTE energy 
bandpass. 

In this paper, we describe the end of J1701's outburst as 
the source transitioned from active accretion to quiescence 
and report on our monitoring of the subsequent cooling of 
the NS in quiescence. 

2. DATA ANALYSIS AND RESULTS 

2.1. RXTE Analysis 

A complete spectral analysis of RXTE da t a from the out- 
burst of J 1701 can be found in iLin et all 1)2009 bP . Here 
we present results fro m the analysis of RX TE Proportional 
Counter Array fPCA; Uahoda et al.|[2006h data from 61 ob- 
servations made during the last phase of the outburst and 
early quiescence (2007 July 17- August 29), with the main 
goal of obtaining flux values in the 0.5-10 keV band. For 
our analysis we only used data from Proportional Counter 
Unit 2 (PCU2). For each of the observations, a single PCU2 
spectrum was extracted from Standard-2 mode data using 
HEASOFT, version 6.8. The spectra were corrected for 
dead time, background was subtracted, and a systematic 
error of 0.6% was added to account for uncertainties in the 
PCA resp onse . The three type I X-ray bursts discussed 
in Section ll.il were removed from the data. A preliminary 
inspection of the light curves shows that the PCU2 count 
rate reached a nearly constant value of ~2 counts s" 1 after 
August 7. This flux is probably due to Galactic background 
emission, since Swift observations made in the last few days 
of the outburst (see Section [22]) indicate much lower fluxes 
than contemporaneous RXTE observations. A spectrum of 
this residual emission was created by combining the spec- 
tra taken between August 8 and August 29. This spectrum 
was then subtracted from the spectra of observations taken 
before August 8; later observations were not considered fur- 
ther. We note that for the August 7 observation the resid- 
ual emission represented ~2/3 of the 3.2-25 keV flux before 
subtraction. 

The sp ectra were fit ted in the 3.2-25 keV range with 
XSPEC ()ArnaudlH996l) . version 12.5.1. For the first ob- 
servations we used a model consisting of a multicolor 
disk blackbody (diskbb in XSPEC), a blackbody (bbody), 
a power law (powerlaw), and a Gaussian emission line 
(gauss). These components were modified by a photoelec- 
tric absorption component (phabs) with the column den- 
sity fixed at Nh = 1-93 x 10 22 cm~ 2 (as determined from 
our fits to the Chand ra and XMM-New ton spectra; see Sec- 
tiondXSl). FollowingEineLaLl (2009b) we fixed the energy 
and width of the Gaussian at 6.5 keV and 0.3 keV, respec- 
tively, and constrained the power-law index to be lower 
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Figure 2. Total unabsorbcd luminosity in the 0.5-10 keV band 
around the end of the outburst. The two lines are best-fit exponen- 
tial decay curves for the three Swift observations, and the first three 
Chandra and XMM-Newton observations. The intersection of these 
curves defines the end time of the outburst, to. 

than or equal to 2.5. In addition to these constraints, we 
also fixed the normalization of the diskbb component at a 
value of 17.0 in XSPEC; this component otherwise becomes 
poorly constrained during the decay, leading to large uncer- 
tainties when extrapolating the spectral model below the 
RXTE energy bandpass. The value of 17.0 is an average of 
values found for the disk normalization in the atoll phase 
of the outburst (when the radius of the disk was approxi- 
mately constant). The diskbb and bbody components were 
statistically no longer needed in spectra taken after approx- 
imately July 29 08:00 and August 4 12:00 UT, respectively, 
and were dropped from the spectral model after those times. 
The best-fit model was extrapolated down to 0.5 keV (well 
below the lower-energy boundary of the PCA) to obtain 
0.5-10 keV unabsorbed fluxes. Such an extrapolation can 
give rise to significant systematic errors in the flux values, 
due to uncertainty in the choice of a proper spectral model. 
This is on top of possible error associated with the subtrac- 
tion of the background Galactic emission. However, the 
good agreement between the RXTE and contemporaneous 
Swift observations (see Figure [2]) indicates that the RXTE 
flux values do not suffer from large systematic errors. The 
derived fluxes for the first Swift observation and a simul- 
taneous RXTE observation differ by only ~6% and agree 
within the uncertainties due to counting statistics. We note 
that since this paper is primarily concerned with the quies- 
cent phase of the source, none of our results are affected by 
the values of the RXTE fluxes. When calculating the lumi- 
nosities for those RXTE observations that included a disk 
blackbody component in their spectra (i.e., observations in 
the first ~13 days in Figure Ht, a disk inclination of 70° 
was assumed (|Lin et al . 2009b). Possible error in these lu- 
minosities due to uncertainty in the assumed inclination 
is probably at most a factor of 2 for the earliest observa- 
tions; less for the later ones, since the contribution of the 
disk blackbody component to the total flux gradually de- 
creases. The RXTE errors plotted in Figure [2] are the same 
fractional errors as those derived for the absorbed 3.5-10 
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Table 1 

Swift, Chandra, and XMM-Newton Observations of XTE J1701-462 



Designation 


Start Date 


Instrument 


ObsID 


Good Exp. Time (ks) 


Sw-1 


2007 


Aug 2 


Swift XRT 


00030383023 


1.82 


Sw-2 


9007 


Aug 4 


Swift XRT 


00030383024 


1.93 


Sw-3 


2007 


Aug 8 


Swift XRT 


00030383026 


2.29 




0007 

ZUU / 


Aug lz 




/ olo 




CXO-2 


2007 


Aug 20 


Chandra ACIS-S 


7514 


8.78 


XMM-1 


2007 


Aug zo 


XMM-Newton EPIC 


0413390101 


20.15 (MOSl) 
19.27 (MOS2) 
9.51(pn) 


XMM-2 


2007 


Sep 28 


XMM-Newton EPIC 


0413390201 


22.02 (MOSl) 
23.12 (MOS2) 
12.36 (pn) 


CXO-3 


2008 


Jan 31 


Chandra ACIS-S 


7515 


19.91 


XMM-3 


2008 


Mar 22 


XMM-Newton EPIC 


0510990101 


30.36 (MOSl) 
31.05 (MOS2) 
21.65 (pn) 


CXO-4 


2008 


Jun 3 


Chandra ACIS-S 


7516 


27.37 


CXO-5 


2008 


Oct 13 


Chandra ACIS-S 


7517 


39.92 


CXO-6 


2009 


Jan 30 


Chandra ACIS-S 


10063 


50.06 


CXO-7a 


2009 


Mar 23 


Chandra ACIS-S 


10064 


13.26 


CXO-7b 


2009 


Mar 23 


Chandra ACIS-S 


10891 


9.02 


CXO-7c 


2009 


Mar 24 


Chandra ACIS-S 


10889 


11.88 


CXO-7d 


2009 


Mar 25 


Chandra ACIS-S 


10890 


19.57 


CXO-8 


2009 


May 23 


Chandra ACIS-S 


10065 


62.31 


CXO-9 


2009 


July 16 


Chandra ACIS-S 


10066 


65.36 


CXO-lOa 


2009 


Oct 13 


Chandra ACIS-S 


10067 


43.56 


CXO-lOb 


2009 


Oct 13 


Chandra ACIS-S 


12006 


25.80 



keV flux in XSPEC, and only take into account uncertainty 
due to counting statistics (in addition to the assumed 0.6% 
systematic error arising from uncertainty in the PCA re- 
sponse). We note here that all errors quoted in this paper 
correspond to la Gaussian (68.3%) confidence. 

2.2. Swift Analysis 

J1701 was observed 3 times w ith the Swift X-Ray Tele- 
scope (XRT; lBurrows et al.ll2005T) in the last few days of the 
outburst (see Table Q]) . All three observations were made 
with the detector in the Photon Counting mode. We ana- 
lyzed the data using HEASOFT (ver. 6.8), and the latest 
calibration files available at the time of the analysis. Start- 
ing with the Level 1 raw event files, the data were processed 
and screened using the xrtpipeline task with the default 
parameter settings, as described in The Swift XRT Data 
Reduction GuideO Additionally, a short period with an in- 
creased background count rate was filtered out in the third 
observation. We extracted source and background spectra 
from the Level 2 screened event files using Xselect. The 
first observation suffered from moderate pile-up and there- 
fore a core of radius 12" was excised from an extraction 
circle of radius 70". The appropriate excision radius was 
determined by comparing the observed point-spread func- 
tion (PSF) to the em pirically determined (pile-up free) one 
(|Moretti et al.l I2005D . For the second and third observa- 
tions, source spectra were extracted from circular regions 
of radii 47" and 35", respectively. Background spectra were 
in all three cases extracted from a source-free circular re- 
gion of radius 200". We used a standard redistribution 
matrix file (RMF) from the calibration database and cre- 
ated observation-specific "empirical" (i.e., adjusted to fit 
the Crab spectrum; see the Data Reduction Guide) ancil- 
lary response files (ARFs) using the xrtmkarf task. The 
ARFs were corrected for counts missed by the finite ex- 

8 Available at http: //swift .gsfc.nasa.gov/docs/swift/analysis/xrt 



traction regions (as were our Chandra and XMM-Newton 
ARFs). 

The spectra were fitted in the 0.2-10 keV band with 
XSPEC (ver. 12.5.1). We binned the spectra from the first 
two observations into groups with a minimum of 25 counts. 
Due to the low number of source counts in the third ob- 
servation (52), that spectrum was left unbinned and fitted 
using the W statistic (a modification o f the C statistic that 
allows for background subtraction; see IWachter et aH[l97fll 
and the XSPEC User's GuideJB). We fixed the absorption 
column at the value of 1.93 x 10 22 cm -2 derived from our 
fits to the Chandra and XMM-Newton spectra. All three 
Swift spectra were adequately fitted with a simple absorbed 
power law. For none of the observations did the addition of 
a thermal component to the spectral model improve the fit. 
The best-fit values for the power-law index are 1.68 ± 0.06 
for the first observation, 1.92 ±0.09 for the second one, and 
2.8 ± 0.4 for the third. 

2.3. Chandra Observations 

All the Chandra observations made after J1701 en- 
tered quiescence (s ee Table [J) used ACIS-S imaging 
(iGarmire et al.l [20031 ) in the Timed Exposure mode, with 
the source located at the nominal aimpoint on the S3 chip. 
The detector was operated either in full-frame (frame time 
of ~3.2 s) or 1/8-subarray mode (frame time of ~0.4 s), us- 
ing the Faint or Very Faint telemetry form at. We analyzed 
the o bservations using the CIAO software (|Fruscione et al.l 
2006), v ersion 4.2 (CALDB ver. 4.2.0), and with ACIS Ex- 
tract (|Broos et al.|[2010h . version 2010-01-07. The data 
were processed following the standard ACIS data prepara- 
tion procedure recommended by the Chandra X-Ray Cen- 
terj__] We checked the observations for possible periods of 

9 Available at http: //heasarc .gsf c .nasa. gov/docs/xanadu/xspec/ 

10 The ACIS Extract software package and User's Guide are avail- 
able at http://www.astro.psu.edu/xray/acis/acis_analysis.html 
swgUiSeQ'Mt_Ep.p<lfcxc .harvard. edu/ciao/guides/acis_dat a. html 
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background flaring, using background light curves from the 
entire SI chip (for the full- frame observations) or the entire 
active area of the S3 chip (for the subarray observations) . 
No such periods were found, except for a short spike in Ob- 
servation ID 10066 which was excluded from the data. We 
performed further analysis with the help of ACIS Extract. 
Source spectra were extracted from polygon-shaped regions 
modeled on the Chandra ACIS PSF using the MARX ray- 
trace simulator. The extraction regions had a PSF enclosed 
energy fraction of ~0.97 (for a photon energy of ~1.5 keV) 
and a radius of ~l'/9. Background spectra were extracted 
from source-free circular annuli centered on the source lo- 
cation, with an inner radius of 1.5 times the radius that 
encloses 99% of the PSF. We chose the outer radius so that 
the background region had at least 4 times the exposure- 
corrected area of the source extraction region, and a mini- 
mum of 100 counts (the latter condition in all cases leading 
to a much higher ratio than 4 between the areas of the 
background and source regions). Response files were cre- 
ated using the mkacisrmf and mkarf tools in CIAO. 

2.4. XMM-Newton Observations 

The three XMM-Newton observations mad e of J 1701 in 
quiescence (see Table [U u sed EPIC imaging (| Turner et al.l 
120011: iStriider et all 1200 lb with the MOS detectors oper- 
ated in the small-window or full-frame mode, and the pn 
detector in the extended-full-frame mode. The thin filter 
was used in all cases. We analyzed the data using the SAS 
software, version 9.0.0, and the latest calibration files avail- 
able at the time of the analysis. Starting with the original 
Observation Data Files, the data were reprocessed using 
the emproc and epproc tasks. Significant portions of the 
exposures had to be excluded due to periods of increased 
particle background. These periods were identified by con- 
structing light curves from the entire area of the detectors, 
using only single events (PATTERN=0) and high energy 
photons (larger than 10 keV for MOS and between 10 and 
12 keV for pn). Intervals with count rates above a cer- 
tain level (0.18-0.24 counts s" 1 for MOS, 0.36-0.40 counts 
s _1 for pn) were then excluded. An additional time cut 
was made to better filter out a large flare at the end of 
the third observation. We used the evselect task to ex- 
tract spectra, after utilizing the eregionanalyse task to 
optimize the circular source extraction region in each case 
(i.e., center the region on the centroid of the counts dis- 
tribution and find the radius that maximizes the ratio of 
source counts to background counts). Extraction radii in 
the range 25"-45" were used, giving PSF enclosed energy 
fractions in the range ~0. 83-0. 87. Background spectra were 
extracted from source-free circular regions with roughly 4 
times the area of the sour ce extraction regions. Follow- 
ing the recommendations in lGuainazzi et al.l (|2010f ). we ex- 
tracted background spectra for the MOS detectors from an- 
other region on the same CCD as the source was on, away 
from source counts; for the pn detector, the background re- 
gion was on an adjacent CCD to the one the source was on, 
at a similar distance to the readout node as the source re- 
gion. The event selection criteria used were PATTERN=0- 
12 and FLAG=#XMMEA.EM for the MOS detectors, and 
PATTERN=0-4 and FLAG=0 for the pn detector. We cre- 
ated both the MOS and pn spectra with a resolution of 30 
eV; this ensured that the grouped spectra used for fitting 
(see Section I2.5.2|) would not oversample the resolution of 



the detectors by more than a factor of 2-3. Response files 
were created using the rmf gen and arf gen tasks. 

2.5. Fitting of Chandra and XMM-Newton Spectra 
2.5.1. Spectral Models 

The Chandra and XMM-Newton spectra were modeled 
with two source components, a NS atmosphere model 
(thermal component) and a simple power-law model (non- 
thermal component), along with an overall photoelectric 
absorption component. For the latter, we used the phabs 
model in XSPEC with the default cross sections and rela- 
tive abundances. 

NS atmosphere spectra have a blackbody-like shape, but 
with an important difference to classic blackbodies being 
a shift in the peak of the emitted radiation to higher pho- 
ton energies and a harder high-energy tail. This is due 
to the strong energy dependence of free-free absorption in 
the NS atmosphere leading to higher energy photons com- 
ing f rom deeper and hotte r layers in the atmosphere (see, 
e.g.. lZavlin fc Pavl ov 2002). Fitting spectra from NS atmo- 
spheres with a classic blackbody model therefore leads to an 
overestimate of the effective temperature, and to an under- 
estimate of the radius o f the emitting region by as much as 
an order of magnitude ( Rutlcdg e et al.lll999f) . We used the 
NS atmosphere model nsatmos, avail able in XSPEC (fo r 
detailed information on the model, see iHeinke et ail feOOe). 
This model has five fitting parameters: the unredshifted ef- 
fective surface temperature of the NS (T e g), the NS mass 
(M ns ), the true NS radius (R ns ), the distance to the NS 
(D) , and an additional normalization parameter represent- 
ing the fraction of the NS surface emitting radiation. For 
given values of M ns and R ns one can, from T ff , calculate the 
more often quoted redshiftcd effective temperature as mea- 
sured by an observer at infinity as = T e ff/(l + z), where 
1 + z = (1 — Rs/ Rns)^ 1 ^ 2 is the usual gravitational redshift 
factor, with Rg — 2GM ns /c 2 being the Schwarzschild ra- 
dius. One can also calculate the observed NS radius at 
infinity as R^ = J? ns (l + z). The canonical values for the 
NS parameters of M ns = 1.4 Mq and R ns — 10 km give a 
redshift factor of 1 + z ~ 1.306. 

The nsatmos model makes several simplifying assump- 
tions, including the following. (1) Magnetic fields are as- 
sumed to have negligible effects on the spectrum from the 
NS atmosphere (i.e., B < 10 8 -10 9 G). Given that J1701 has 
shown behavior typical of low-magnetic-field NS-LMXBs, 
with no X-ray pulsations having been detected, this should 
be a reasonably good assumption. (2) The NS atmo- 
sphere is assumed to be composed of pure hydrogen. Since 
heavier elements are expecte d to settle out of the atmo- 
spher e on a timescale of ~10 s (|Romanill9 87: Bildstcn et al.l 
1992), this should be true if the accretion rate is below 
~10 -13 M Q yr" 1 (|Brown et al.lll998D . If the non-thermal 
component seen in the spectra from J1701 is due to a resid- 
ual accretion flow which is radiating efficiently, then the 
accretion r ate is likely ~10~ 14 -10~ 12 Mq yr _1 for our ob- 
servations. IHeinke et alj (|2006T ) note that atmospheres with 
iron or with solar abundances would be easily identifiable 
by their different spectral shapes, whereas small departures 
from a pure hydrogen atmosphere may go unnoticed and af- 
fect results. It is therefore conceivable that our results may 
be affected by some contamination of the hydrogen atmo- 
sphere, although the fact that our spectra are in general well 
fitted with nsatmos (with an implied emission area consis- 
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Table 2 

Selected Spectral Fit Parameters 



Fit 


D (kpc) 


R nB (km) 


M ns (M ) 


Nn (10 22 cm" 2 ) 


a a 


a 'b 


1 


(8.8) 


(10) 


(1.4) 


1.93 ± 0.02 


1.93 ±0.20 


1.35 ±0.08 


2 


(7.5) 


(10) 


(1.4) 


1.98 ±0.03 


2.04 ±0.19 


1.37 ±0.08 


3 


(10.1) 


(10) 


(1.4) 


1.88 ±0.02 


1.83 ±0.20 


1.34 ±0.08 


1 


7.3 ±1.4 


(10) 


(1.4) 


1.99 ±0.07 


2.06 ±0.23 


1.37 ±0.08 


5 


(8.8) 


11.7 ±2.0 


(1.4) 


1.98 ±0.06 


2.02 ± 0.22 


1.36 ±0.08 


6 


(8.8) 


(10) 


1.6 ±0.4 


1.95 ± 0.06 


1.98 ±0.21 


1.36 ±0.08 


7 


(8.8) 


11.6 ±2.6 


1.5 ±0.6 


1.98 ±0.07 


2.03 ±0.24 


1.36 ±0.08 


8 


(8.8) 


(10) 


(1.4) 


1.95 ±0.02 


(1.3) 


1.34 ±0.08 


9 


(8.8) 


(10) 


(1.4) 


1.92 ± 0.02 


(2.5) 


1.35 ±0.08 



Note. — All Chandra and XMM-Newton observations were fitted simultaneously. All 
the fits have a reduced x 2 of 1.11-1.13 for 515—518 dof. Numbers in parentheses indicate 
parameters fixed during fitting. Errors quoted are at the la Gaussian (68.3%) confidence 
level. — a Combined (tied) power-law index for all observations except the third XMM- 
Newton observation (XMM-3). b Power-law index for XMM-3. 



tent with what is expected from a NS) suggests that this is 
not a serious problem. (3) Finally, it is assumed that the 
hydrogen is completely ionized, and Comptonization is ig- 
nored. This limits the validity of the co de to a temperatur e 
range of 3 x 10 5 K < T eS < 3 x 10 6 K (IHeinke et al.ll2006l) . 
corresponding to 20 eV < kT^ < 200 eV for typical values 
of the NS parameters. The derived effective temperatures 
for J 1701 fall well within this range. 

In addition to nsatmos, two other spectral models for NS 
atmospheres with negligible magnetic fields are available 
in XSPEC, nsa and nsagrav. The nsatmos and nsagrav 
models allow for variations in surface gravity correspond- 
ing to different values of the NS radius and mass. In con- 
trast, nsa uses a single fixed surface gravity value. Com- 
parisons of different NS atmosphere models have shown 
that taking variations in surface gravity into account can 
be im portant when allow i ng the NS mass and/o r radius to 
vary (Hcinkc eTaIll2006t I Webb fc Barret! 120071) : nsatmos 
and nsagrav are therefore to be preferred. We did a cur- 
sory comparison of our spectral fitting results when using 
nsatmos and nsagrav, and found them to be largely equiv- 
alent. In what follows we only report results obtained using 
nsatmos. 

We model the non-thermal component with a simple 
power law (using the pegpwrlw model in XSPEC, whose 
normalization is pegged to equal the energy flux of the 
model for a chosen energy range). We note that this is 
not a physically motivated choice (the origin of this com- 
ponent is uncertain), but is simply based on the fact that 
non-thermal components in quiescent NS-LMXB spectra 
have often been successfully modeled by power laws in the 
past. It is, however, very possible that our results may be 
somewhat affected by the power-law model not properly 
representing the true shape of the non-thermal component. 
To get some indication of how our results might be affected 
by the choice of a model for the non-thermal component, 
we also performed a fit where the power-law model was re- 
placed with the simpl model in XSPEC (see discussion in 
Section [2HT]). 

2.5.2. Fitting Results 

We fitted the Chandra and XMM-Newton spectra in the 
0.5-10 keV band with XSPEC (ver. 12.5.1), after binning 
the spectra into groups with a minimum signal-to-noise ra- 
tio of 4 using the ACIS Extract tool ae_group_spectrum 
(resulting in an average number of counts per group of 



~25-30 for the Chandra spectra and ~35-40 for the XMM- 
Newton spectra). The spectra from all 13 Chandra and 
XMM-Newton observations were fitted simultaneously. The 
seventh Chandra observation (CXO-7) consisted of four sep- 
arate exposures taken over a period of approximately 3 
days (see Table [I}. In the main spectral fit, used to de- 
rive temperatures and fluxes, we tied all the parameters for 
these four spectra, thereby effectively treating the expo- 
sures as a single observation. However, we also performed 
a fit without the parameters being tied (see Section 12.6. ip . 
The tenth and most recent Chandra observation consisted 
of two separate exposures taken over a ~27 hr period; this 
observation was treated in the same manner as CXO-7. For 
each XMM-Newton observation, all fitting parameters were 
tied between the spectra from the three EPIC detectors 
(MOS1, MOS2, and pn). Several parameters were tied be- 
tween all 13 observations: the absorption column, the NS 
mass, radius, and distance, and the fraction of the NS sur- 
face emitting. Additionally, due to the limited statistics 
of our spectra and the fact that we are mainly interested 
in constraining the values for the effective temperature (as 
well as the flux contribution from the non-thermal compo- 
nent), most of the nsatmos parameters were usually fixed 
at a certain value and not allowed to vary during the fit- 
ting. (1) The value of the distance parameter was fixed at 
the best-estimate value of 8.8 kpc (see § ll.lj) in the main 
fit, although we also performed fits where the distance was 
allowed to vary or was fixed at 7.5 kpc or 10.1 kpc. (2) 
Likewise, the NS mass and radius were fixed at the canon- 
ical values of 1.4 M Q and 10 km in the main fit, although 
we also experimented with allowing them to vary (see be- 
low). (3) Finally, since we assume that the entire surface 
of the NS is emitting thermal radiation during quiescence, 
the fraction of the NS surface emitting was always fixed at 
1. 

Although we have no reason to believe that the shape of 
the non-thermal component should necessarily be the same 
in all the observations (as mentioned above, its origin is 
highly uncertain), we do not have enough counts to allow 
the power-law index to vary freely for each observation. 
Doing so leads to very poorly constrained values for the in- 
dex in most cases and a very unstable fit. An exception is 
the third XMM-Newton observation (XMM-3), which has 
a much larger non-thermal component than the other ob- 
servations. In the main fit, we therefore tie the power-law 
index between all the observations except XMM-3. Ty- 
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Table 3 

Derived Temperatures and Fluxes from Chandra and XMM-Newton Observations 



\J DSGrVeltlOIl 


f v-„ a 

l l() 

( H a vs I 


(eV) 


-^bol 

(10 — ^ 3 erg s — ^ cm - ^) 


F , d 

(10 Grg s ^ cui - ^) 


Mot 
fin 33 ere- s - M 


CXO-1 


2.77 


164.2 ±3.6 


17.2 ± 1.5 


1.5 ±1.1 


15.2 ± 1.7 




10.63 


1 £n k _l_ o £ 
159.0 ± 2.5 


15.4 ± 1.0 


1.3 ± 0.7 


13.4 ± 1.2 


XMM-1 


16.06 


156.8 ± 1.3 


14.3 ±0.5 


0.5 ±0.3 


11.8 ± 0.5 


XMM-2 


49.31 


150.0 ±1.2 


12.0 ±0.4 


0.8 ±0.3 


10.1 ± 0.5 


CXO-3 


174.15 


129.1 ±4.7 


6.6 ± 1.0 


4.8 ±0.9 


q .,+0.6 


XMM-3 


225.54 


159.3 ±2.0 


15.3 ±0.8 


15.1 ±0.6 


26.1 ± 0.5 


CXO-4 


298.12 


136.0 ±2.0 


8.1 ±0.5 


1.4 ±0.3 


7.3 ±0.5 


CXO-5 


430.89 


126.3 ±3.1 


6.0 ±0.6 


3.1 ±0.5 


7 O+0.4 
'■ Z -0.5 


CXO-6 


539.90 


125.4 ±1.5 


5.8 ±0.3 


0.5 ±0.2 


4.7 ±0.3 


CXO-7 f 


592.50 


129.6 ±2.2 


6.7 ±0.5 


2.0 ±0.4 




CXO-8 


652.44 


124.0 ±2.2 


5.6 ±0.4 


1.8 ±0.3 


t- 7+O.3 


CXO-9 


705.20 


123.9 ±2.0 


5.6 ±0.4 


1.5 ±0.3 


5 4+0.3 


CXO-10 g 


795.45 


124.1 ± 1.7 


5.6 ±0.3 


1.0 ±0.2 


5.0 ±0.3 


Note. — Values were 


obtained from 


a simultaneous spectral fit to all observations, with the NS mass, 



radius, and distance fixed at 1.4 Mq, 10 km, and 8.8 kpc (fit 1 in Table [2j. Errors quoted are at the 
lcr Gaussian (68.3%) confidence level. — a Time of mid-observation; to is MJD 54322.13. b Effective 
surface temperature of the NS as seen by an observer at infinity. c Unabsorbed bolometric flux of the 
thermal (nsatmos) component. d Unabsorbed 0.5-10 keV flux of the non-thermal (power-law) component. 
e Unabsorbed total luminosity in the 0.5—10 keV band. f Values derived with parameters for the CXO- 
7[a-d] spectra tied. g Values derived with parameters for the CXO-10[a-b] spectra tied. 

acceptable fit, and is included in all observations for con- 
sistency. 

In Table [5J we show the results of spectral fits with differ- 
ent combinations of fixed and free NS parameters. We do 
not show the parameters that were allowed to vary freely 
for each individual observation (i.e., the effective tempera- 
tures and power-law normalizations); those would add 26 
additional numbers for each fit to the table. Fit 1 in Ta- 
ble [2] is the main fit that was used to derive temperatures 
and fluxes used in most of the subsequent analysis; these 
are shown in Table [3] Spectra from this fit for two of 
the observations (those with the highest and lowest total 
fluxes) are shown in Figure [H The other fits in Table [2] 
explore other values for the NS parameters and tied power- 
law index. Allowing the distance to vary (fit 4) for fixed 
canonical values of the NS mass and radius yields a best-fit 
distance of 7.3 ± 1.4 kpc, con sistent with the distance esti- 
mate from ILin et al.l (|2009aD . Allowing instead the radius 
to vary (fit 5) gi ves R ns = 11.7 ± 2.0 km. We note that 
ILin et all (|2009aD independently derive an asymptotic (ap- 
parent blackbody) radius of 8 ± 1 km from the three type I 
bursts and, correcting for the effects of redshift and spectral 
hardening (assuming a distance of 8.8 kpc, a NS mass of 
1.4 Mq, and a hardening factor of 1.4), quote an actual NS 
radius of ~13 km. This value must, however, be regarded 
as highly uncertain, mainly due to the large uncertainty in 
the appropriate value for the hardening factor. Allowing 
the mass to vary (fit 6) gives M ns = 1.6 ± 0.4 Mq. Allow- 
ing both the radius and mass to vary (fit 7) gives a radius 
of 11.6 ± 2.6 km and a mass of 1.5 ± 0.6 Mq. The derived 
absorption column and power-law indices are largely unaf- 
fected by the values of the other parameters. The values of 
the absorption column and the power-law index for XMM-3 
are tightly constrained to be Ah — (1.9-2.0) x 10 22 cm~ 2 
and a' ~ 1.3-1.4. This value for Ah is in good agreement 
with the value of ~2.0 x 10 22 cm -2 deri ved from RXTE a nd 
Swift observations during the outburst (|Lin e t al. 200930). 
The combined (in some sense averaged) power-law index 
for the other observations (i.e., all the quiescent observa- 
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Figure 3. Fits to the pn spectrum for the third XMM-Newton obser- 
vation (upper) and the spectrum for the sixth Chandra observation 
(lower). Also shown (dotted lines) are the contributions from the 
two spectral components: the thermal component, which dominates 
at lower energies, and the non-thermal component, which dominates 
at higher energies. The upper spectrum has the largest non-thermal 
component of all the observations, whereas the lower spectrum has a 
non-thermal component among the smallest seen. 

ing the XMM-3 power-law index to the other observations 
would lead to it completely dominating the fit for that pa- 
rameter. To gauge the effect of this tying on the derived 
temperatures and fluxes we also performed fits where the 
index was fixed at 1.3 and 2.5 for all observations except 
XMM-3 (see discussion in Scction r2.6.1[) . The only param- 
eters allowed to vary independently for each observation in 
the main fit were the effective NS surface temperature and 
the power-law normalization. We note that the values ob- 
tained for the power-law indices when allowing them to vary 
freely for each observation are nearly all consistent (within 
their large error bars) with the single value obtained when 
tying the index as described above. We also note that the 
power-law component is clearly required for nearly all the 
individual observations (and not just XMM-3) to get an 
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Figure 4. Total unabsorbed luminosity in the 0.5-10 keV band (top 
panel), redshifted effective NS surface temperature (middle panel), 
and unabsorbed power-law flux in the 0.5—10 keV band (bottom panel) 
during quiescence. The solid curve in the temperature panel is the 
best-fit exponential decay cooling curve (with the sixth and seventh 
data points excluded from the fit), and the dashed line represents the 
best-fit constant offset to the decay. 



tions except XMM-3) is a less meaningful quantity, but 
seems to be higher than the XMM-3 index, having a value 
of a ~ 1.9±0.3. The main fit has a reduced x 2 value, x 2 > 
of 1.12 for 517 degrees of freedom (dof), corresponding to 
a x 2 probability, P x , of 0.035. This is only a marginally 
adequate fit; possible reasons for this will be mentioned in 
Section POTTI The other fits in Table [2] all have xl in the 
range 1.11-1.13 for 515-518 dof. 

From our main spectral fit we derive both unabsorbed and 
absorbed total fluxes, as well as unabsorbed fluxes for each 
of the two spectral components. We use the same fractional 
errors for the unabsorbed total fluxes as for the absorbed 
ones, thereby ignoring some error arising from uncertainty 
in the effects of the absorption on the flux values. However, 
given how tight our constraints on the absorption column 
are, this should be a relatively small effect. Since the tem- 
perature is the only free parameter in the nsatmos model, 
we can simply propagate the error in the temperature to 
get the error in the unabsorbed thermal flux. The value 
of the normalization parameter of the pegpwrlw model is 
the unabsorbed flux of that component (in a chosen energy 



range), and the error in that parameter therefore directly 
gives the error in the unabsorbed power-law flux. 

2.6. Cooling Curves 

Figure [2] shows the transition from the final stage of out- 
burst to quiescence. Plotted is the total unabsorbed lumi- 
nosity in the 0.5-10 keV band for the 37 RXTE observations 
made in the period 2007 July 17- August 7, and the three 
Swift observations discussed above, as well as the first three 
Chandra and XMM-Newton observations. The luminosity 
decreased by a factor of ^2000 in the final ~13 days of the 
outburst before starting a much slower decay. This period 
of low-level and slowly changing (compared to the outburst 
phase) emission, taking place after the steep drop in lu- 
minosity, is what we refer to as the quiescent phase (see 
also the top panel in Figure |4|) . Low-level accretion may 
be occurring during quiescence, but this current phase is 
clearly distinct from the much more luminous and variable 
outburst phase, during which accretion took place at much 
higher rates (and which we also refer to as the period of 
"active" accretion). The end of the outburst is tightly con- 
strained to have occurred sometime in the ~4.3 day interval 
between the final Swift observation and the first Chandra 
observation. To get a more precise estimate for the end 
time of the outburst, here denoted by to, we fit simple ex- 
ponential decay curves through the three Swift data points 
and the three Chandra and XMM-Newton points in Fig- 
ure [2] From the intersection of those two curves we define 
t as MJD 54322.13 (2007 August 10 03:06 UT), i.e., -2.8 
days before the first Chandra observation. 

Table [3] lists temperatures and fluxes derived from the 
main fit to the Chandra and XMM-Newton spectra dis- 
cussed in Section [2321 Figure |4] shows a plot using results 
from this fit. The top two panels show the total unabsorbed 

0. 5-10 keV luminosity and the inferred effective NS surface 
temperature (as observed at infinity). The first five data 
points, taken in the first ~175 days of quiescence, show 
a fast drop in temperature. However, the sixth data point 
(XMM-3; at ~226 days) shows a large increase in both tem- 
perature and luminosity, and the following Chandra obser- 
vation (CXO-4) also has a higher inferred temperature than 
before the increase. This is inconsistent with the mono- 
tonic decrease in temperature expected for a cooling NS 
crust. The last six Chandra observations all have tempera- 
tures similar to or slightly lower than the one immediately 
preceding XMM-3 (i.e., CXO-3). We assume that those 
are unaffected by whatever caused the "flare-like" behav- 
ior in the sixth and seventh observations, and when fitting 
cooling models to the data we exclude both XMM-3 and 
CXO-4 but include the subsequent observations (although 
some fits excluding only XMM-3 were also made; see be- 
low). We defer further discussion of the flare to the end of 
this section and Section l3~3l 

We will now describe our fitting of the derived temper- 
atures with cooling curve models. All the fits were per- 
formed with Sherva, CIAO's modeling and fitting package 
(jFreeman eEaD[2001); errors were estimated with the confi- 
dence methodQ We first fitted our temperature data with 
an exponential decay cooling curve plus a constant offset, 

1. e., a function of the form T^(t) = V exp[-(t-t )/T]+T eq , 
with to kept fixed at the value mentioned above. Shifts in 



12 See documentation a t 
http : //exc . harvard . edu/sherpa/ 



the Sherpa Web site: 
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Table 4 

Best-fit Cooling Curve Parameters for Temperature Data 



Spectral Fit Parameters 



Exponential Decay Fit 



Broken Power-law Fit 



D 

(kpc) 


Has 

(km) 


(Mq) 


a a 


(days) 


(eV) 


kT' a 
(eV) 


AC 


71* 


72 g 


/ h 

*b 

(days) 




(8.8) 


(10) 


(1.4) 


1.93 




125.0 ±0.9 


37.4 ± 1.8 


0.84 


0.027 ±0.013 


0.070 ± 0.004 


QQ+23 


1.08 


(7.5) 


(10) 


(1.4) 


1.93 


118±?5 


118.4 ±0.9 


35.2 ± 1.7 


0.76 


0.025 ±0.013 


0.070 ± 0.004 


32+f 


1.04 


(10.1) 


(10) 


(1.4) 


1.93 




131.0 ±0.9 


39.5 ± 1.8 


0.97 


0.030 ±0.013 


0.070 ± 0.004 


34+^ 


1.15 


(8.8) 


(11.6) 


(1.5) 


1.93 


117+27 
111 -19 


120.3 ±0.9 


35.9 ± 1.7 


0.76 


0.025 ±0.013 


0.070 ± 0.004 


32+f 


1.04 


(8.8) 


(10) 


(1.4) 


(1.3) 


1 QO + 27 

1,3z -20 


128.2 ±0.7 


34.1 ± 1.3 


2.89 


0.030 ±0.011 


0.063 ± 0.003 




2.46 


(8.8) 


(10) 


(1.4) 


(2.5) 




120.8 ±0.9 


41.5 ±2.1 


1.75 


0.023 ±0.017 


0.079 ± 0.005 


30tf 


2.61 



Note. — In each case, the data points used for the cooling curve fit were derived from a simultaneous spectral fit to all 
observations. Selected spectral fit parameters are shown in the first four columns; numbers in parentheses indicate parameters 
fixed during the spectral fitting. Errors quoted are at the la Gaussian (68.3%) confidence level. — a Combined power-law 
index for all observations except the third XMM-Newton observation (XMM-3). b Best- fit e- folding time of the decay. c Best-fit 
constant offset to the decay. d Best-fit normalization coefficient of the decay. e Reduced x 2 f° r the fit, which had 8 dof in each 
case. f Best-fit pre-break power-law slope. g Best-fit post-break power-law slope. h Best-fit break time between power laws. 
1 Reduced x 2 f° r the fit, which had 7 dof in each case. 
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Figure 5. Effective NS surface temperature during quiescence, with 
best-fit cooling curves shown. The solid and short-dashed curves are 
broken power laws fitted to data points 1-5 plus 8-13, and data points 
1-5, respectively. The gray long-dashed curve is the best-fit exponen- 
tial decay curve with a constant offset (also shown in Fig. [4j . 

the value of to do not affect derived values for r or T eq . The 
flare observations XMM-3 and CXO-4 were excluded from 
the fitting. We performed the temperature fit for data from 
the main spectral fit (1 in Table [2]), and also for spectral 
parameter values corresponding to five other fits (2, 3, 7, 8, 
and 9), to gauge the effects on the cooling fit parameters. 
The derived parameter values are shown in Table |H The 
main fit cooling curve is shown in Figure HI along with the 
best-fit constant offset (dashed line). The best- fit e- folding 
time is r = 117±H days with an offset of T oq = 125.0 ± 0.9 
eV. For the other values of the NS parameters (mass, ra- 
dius, and distance), the temperature values are systemat- 
ically shifted by typically 5-10 eV, but the derived decay 
timescale is not affected to a significant extent. The effects 
of changing the value of the tied power-law index will be dis- 
cussed in Section [2.6.1l Including CXO-4 in the fit (but still 
excluding XMM-3) gives a longer timescale of r — 187^3g 
days; the equilibrium temperature is not significantly af- 
fected. 



As will be discussed in Scction [3.21 a more physically mo- 
tivated cooling curve model is a broken power law leveling 
off to a constant at late times. We therefore also fitted 
a broken power-law model, excluding XMM-3 and CXO-4 
as before, to temperature data corresponding to the same 
six spectral fits as before. The derived break times and 
power- law slopes are shown in Table 0] The best-fit broken 
power-law curve to data from the main spectral fit is shown 
in Figure [5] (solid curve). The data indicate that a break in 
the model is needed; a simple power law does not provide 
an adequate fit (xl = 2.45 for 9 dof, compared to \t = 1-08 
for 7 dof in the broken power- law case) . The best-fit break 
time is t\ } = 33^^ days post-outburst; the best-fit slopes 
are 71 = 0.027 ± 0.013 (pre-break) and 72 = 0.070 ± 0.004 
(post-break) . These values do not change significantly when 
the other values of the spectral parameters are used (see 
Table |4|). The exponential decay and broken power-law fits 
have xt — 1 when the combined power-law index is free 
to vary in the spectral fit; the values for the exponential 
fits tend to be slightly lower. For the spectral fits with 
the power-law index fixed, the xt values are considerably 
higher (~2-3); this is in large part due to the smaller error 
bars for the temperatures resulting from the fixing of the 
index. Including CXO-4 does not significantly affect the 
broken power-law fit. 

Given that the broken power-law cooling curve is ex- 
pected to level off at some point, and that from CXO-3 
onward the observations only show small changes in tem- 
perature (not counting the flare observations), we also fit- 
ted a broken power law to only the first five data points 
from the main spectral fit (short-dashed curve in Figure [SJ. 
This gave the same pre-break slope as before, a break time 
t h = 41±g 31 days, and a post-break slope 72 = 0.119±0.030. 
The upper limit to the break time is in this case only con- 
strained by the time of the fifth observation. The time of 
the break is in general poorly constrained by our fits. In all 
cases, the % 2 surface in parameter space, although locally 
quadratic around the minimum, shows significant devia- 
tions from a parabolic shape toward higher break times, 
in some cases before a difference in the value of \ 2 cor ~ 
responding to a la error is reached. Any upper limits to 
ib given should therefore be regarded with caution. Over- 
all, we conclude that the time of the break could be as 
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Table 5 

Best-fit Cooling Curve Parameters for Bolometric Thermal Luminosity Data 



Spectral Fit Parameters 



Exponential Decay Fit 



Broken Power-law Fit 



D 

(kpc) 


Rns 

(km) 


M ns 
(Mq) 


a a 


r b 
(days) 


T c 

(10 33 erg s" 1 ) 


(10 33 erg s" 1 ) 


AC 

A/7 


71* 


72 s 


(days) 




(8.8) 


(10) 


(1.4) 


1.93 


S8+H 


5.35 ±0.14 


9.83 ±0.61 


0.80 


0.108 ±0.053 


0.282 ± 0.017 


oo + 21 
oo -10 


1.14 


(7.5) 


(10) 


(1.4) 


1.93 


S8+\l 


5.94 ±0.16 


10.80 ±0.68 


0.71 


0.100 ±0.053 


0.281 ± 0.017 


32+1 9 


1.12 


(10.1) 


(10) 


(1.4) 


1.93 


87ti! 


4.90 ±0.12 


9.10 ±0.56 


0.92 


0.119 ±0.053 


0.283 ± 0.016 


34+23 


1.18 


(8.8) 


(11.6) 


(1.5) 


1.93 


S8+\l 


5.87 ±0.16 


10.71 ±0.67 


0.72 


0.101 ±0.053 


0.281 ± 0.017 


32+1 9 


1.13 


(8.8) 


(10) 


(1.4) 


(1.3) 




5.92 ±0.11 


9.15 ±0.50 


2.89 


0.119 ±0.044 


0.255 ± 0.013 


35+?? 


2.37 


(8.8) 


(10) 


(1.4) 


(2.5) 




4.65 ±0.14 


10.73 ±0.69 


1.92 


0.091 ±0.067 


0.316 ± 0.020 


28+f 


3.76 



Note. — All notes for Table l4l apply here as well. 

early as c±20 days and as late as ~150 days post-outburst. 
Since a break time after XMM-2 (which is at ~49 days 
post-outburst) is clearly also allowed by the data, we fit- 
ted a simple power law to the first four data points to 
estimate what the pre-break slope would be in that case. 
This gave 71 = 0.035 ± 0.007, slightly higher than for fits 
with a break before XMM-2. We note that, in contrast to 
the exponential decay fits, the results from the power-law 
fits depend (sensitively) on the value used for to. Plac- 
ing to 2.5 days later (i.e., ~0.3 days before the first ob- 
servation in quiescence) and fitting the entire data set (ex- 
cept XMM-3 and CXO-4) gives 71 = 0.012 ± 0.006 and 
ti, = 26]tg days. However, results from fits with to very 
close to the first data point should be regarded with cau- 
tion. Due to the divergence of the power-law model at to 
(since T°g(t) oc (t — to) -71 ), the power-law behavior of the 
cooling curve cannot extend all the way back to to', the 
curve has to flatten out from a power law before that. We 
also performed a fit with to placed 13 days earlier, i.e., 
around the time when the decay of the outburst light curve 
steepened dramatically (see Figure 0). Although the source 
was still in outburst at that time, it may be appropriate to 
consider the effective start time of the cooling as the onset 
of the fast decrease in accretion rate implied by the steep 
decay. In this case, there is no indication of a break in the 
cooling curve; the entire curve is well fitted with a single 
power law with slope 0.073 ± 0.003, even when including 
CXO-4. 

In addition to the fits to the temperature data discussed 
above, we also made analogous fits to the unabsorbed bolo- 
metric luminosity of the thermal component. This was 
done to permit possible comparison with results from other 
sources, and because such fits are of some interest to theo- 
rists. The results are given in Table [5] and will not be dis- 
cussed further in this paper, apart from noting that the ex- 
trapolation of the nsatmos model outside the energy range 
covered by our data gives rise to some systematic uncer- 
tainty. Our observed 0.5-10 keV thermal flux represents 
~80%-85% of the inferred bolometric thermal flux (practi- 
cally all of which is contained in the 0.01-10 keV band). We 
therefore expect the induced systematic error in the bolo- 
metric thermal flux due this extrapolation to be at most 
-10%. 

The non-thermal flux has varied irregularly throughout 
the quiescent phase (see Table [3] and the bottom panel 
of Figure HJ). As mentioned before, XMM-3 has by far 
the largest non-thermal flux, ~1.5 x 10~ 12 erg s _1 cm~ 2 
in the 0.5-10 keV band, corresponding to an unabsorbed 
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Figure 6. Fractional contribution of the power-law component to 
the total unabsorbed flux in the 0.5-10 keV band as a function of 
total unabsorbed luminosity in the same energy range. The numbers 
indicate the time sequence of the observations. 



luminosity of ~1.4 x 10 erg s for a distance of 



kpc. 

10~ 14 erg s _1 



Values for the other observations range from 

13 „„„ o-l „™-2 



±5 x 



5 x 10 13 erg s 1 cm 2 . A plot of 



the fractional contribution of the non-thermal flux to the 
total unabsorbed 0.5-10 keV flux as a function of the to- 
tal unabsorbed 0.5-10 keV luminosity is shown in Figure|Hl 
We note that changes in the NS parameters in the spectral 
fit yield systematic shifts of ~(l-4) x 10~ 14 erg s _1 cm~ 2 in 
the derived power-law fluxes. The fractional changes are in 
the range ~7%-40% (except for XMM-3 where the changes 
are ~l%-2%). The effects of changing the value of the tied 
power-law index will be discussed in Section 12.6.11 



2.6.1. Further Tests of Cooling Results 

Since we use data from both Chandra and XMM-Newton, 
cross-calibration error between the instruments can possi- 
bly affect our results, especially given how important the 
second XMM-Newton observation is in determining the 
timescale of the decay. The cross-calibration error between 
Chandra ACIS and XMM-Newton EPIC fluxes is in gen- 
eral expecte d to be 5,10% (H. Marshall 2008, priva te com- 
munication; iSnowdenl l2002t iStuhlinger et aLlfeOOSl) . Given 
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the fourth power temperature dependence of the bolomet- 
ric thermal flux, the relative cross-calibration error in the 
temperatures is likely at most a few percent. To estimate 
the effect on our results of a systematic shift in the XMM- 
Newton temperatures with respect to the Chandra ones, 
we increased/decreased the XMM-Newton temperatures by 
±3% and repeated the fits. For the exponential decay 
fit, this yielded e- folding times of r = 149 1 24 days and 

r = 87^ 13 days, respectively, and no significant change in 
the equilibrium temperature. In both cases the quality of 
the fit was worse than without any shifting. The shifts had 
a greater impact on the broken power-law fits. A +3% shift 
gave 71 = 0.007 ± 0.013, 72 = 0.082 ± 0.004, and t h = 30^ 
days, with a decrease in the quality of the fit. After a - 
3% shift there is no longer need for a break in the curve; 
the entire cooling curve is well fitted with a single power 
law with slope 0.054 ± 0.002 (still excluding XMM-3 and 
CXO-4 from the fit, although CXO-4 is in reasonably good 
agreement with the best-fit curve, and including it does not 
change the estimate of the slope) . 

To get an indication of how sensitive our results are to 
our choice of a model for the non-thermal component, we 
repeated our main spectral fit with the power-law model 
replaced by the simpl model in XSPEC. This is an empiri- 
cal convolution model for Comptonization, wh ich converts a 
fracti on of input seed photons to a power law (|Steiner et al.l 
2009). The model has only two free parameters, the power- 
law index and the fraction of scattered photons, and can 
be used with any spectrum of seed photons. Compared to 
our previous results, we see small temperature shifts in the 
range —0.8 to 1.5 eV for all but three of the observations; for 
XMM-3, CXO-3, and CXO-5 (the three observations with 
the largest non-thermal components) we get larger shifts 
of 14.7, 5.5, and 3.5 eV, respectively. The values of the 
absorption column and power-law indices are similar to be- 
fore, as are the xt values. The behavior of the non-thermal 
component also seems to be similar, although comparison 
is complicated by the fact that when using simpl the ther- 
mal and non-thermal components cannot be disentangled in 
the same manner as when a regular power law is used. The 
derived e-folding time for these new temperature values is 
r = 128^21 days with a best-fit constant offset of 126.0±0.8 
eV; these values are not significantly different from the pre- 
vious ones. The quality of the exponential cooling curve 
fit is significantly worse than before, with xt — 1-74 for 8 
dof. In the case of the broken power-law fit, using simpl 
instead of a power-law model does not significantly change 
the derived parameter values, but does adversely affect the 
quality of the fit. 

As mentioned in Section 12.5.21 the seventh Chandra ob- 
servation consisted of four separate exposures taken over 
approximately 3 days, and the tenth Chandra observation 
consisted of two exposures taken over a period of ~27 hr. 
To look for possible variability between exposures within 
these two observations, we performed a spectral fit to all 
the Chandra and XMM-Newton observations, identical to 
the main fit described in Section 12.5.21 apart from the fact 
that the temperature and power-law normalization param- 
eters were not tied between the four CXO-7 spectra and 
the two CXO-10 spectra. Untying these parameters within 
CXO-7 and CXO-10 led to a considerable improvement of 
the overall fit, from xt = 1-12 for 517 dof to xl = 1-07 
for 509 dof (P x = 0.14). Figure [7] shows a plot of the 
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Figure 7. Total unabsorbed luminosity in the 0.5-10 keV band 
(top/first panel), total absorbed flux in the 0.5—10 keV band (sec- 
ond panel), effective NS surface temperature (third panel), and unab- 
sorbed power-law flux in the 0.5-10 keV band (bottom/fourth panel) 
for the four sub-exposures in the seventh Chandra observation. The 
horizontal error bars indicate the durations of the exposures. 

temperature and various fluxes for the four exposures in 
CXO-7. Although not shown in the plots, we note that the 
count rate exhibits a decrease similar to that seen for the 
absorbed flux; the net count rates in the first and last ex- 
posures (in the 0.5-10 keV band) are (1.91 ± 0.12) x 10~ 2 
and (1.57±0.09) x 10~ 2 counts s -1 , respectively. The plot 
indicates a possible decrease in the power-law flux over 
the 3 day period. It is not implausible that this observed 
decrease is real, especially if the non-thermal component 
arises from accretion. However, the fact that the behav- 
ior of the power-law flux seems to be anti-correlated with 
the behavior of the temperature (which shows a slight but, 
given the errors, non-significant increase) may point to a 
limitation in our ability to separate the contributions from 
the thermal and non-thermal components to the total flux. 
Looking at the CXO-10 exposures further indicates that 
this may be the case. The net count rates observed are 
(1.26 ± 0.06) x 10~ 2 and (0.95 ± 0.06) x 10~ 2 counts s" 1 
in the first and second exposures. A short-term light curve 
suggests an overall decrease in net count rate during the 
first exposure; a one-sided Kolmogorov-Smirnov (K-S) test 
comparing photon arrival times to a uniform count rate 
model rejects the null hypothesis of a constant light curve 
at a 97.7% level. However, in this case the spectral fits show 
no change in the power-law flux, but a decrease in temper- 
ature from 127.2 ± 1.9 to 119.5 ± 2.6 eV. It is clearly not 
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plausible for the surface temperature to change in such a 
short time; any real change in flux is likely due to the non- 
thermal component (unless, perhaps, some of the thermal 
flux arises from accretion, see discussion in Section I3.3j) . 
Indeed, constraining the two temperatures to be the same 
and only allowing the power-law normalization to vary in- 
dependently results in a perfectly adequate fit to the two 
spectra; this is also the case for the four spectra in CXO- 
7. It is therefore likely that any observed differences in 
temperatures within CXO-7 and CXO-10 are simply an ar- 
tifact of the spectral fitting, and it is possible that simi- 
lar effects are in general affecting our derived temperatures 
and fluxes to some extent. We note that for the overall fit 
where the temperatures were tied within CXO-7 and CXO- 
10, but the power-law normalizations allowed to vary, the 
derived temperatures are 129.0 ± 2.2 eV for CXO-7 and 
124.3 ± 1.6 eV for CXO-10; these values are not signifi- 
cantly different from the ones in the main spectral fit. In 
this case, the 0.5-10 keV unabsorbed power-law flux for the 
first and last exposures in CXO-7 is (3.2 ±0.5) x 10~ 13 and 
(1.5 ± 0.3) x 10~ 13 erg s" 1 cm" 2 ; for the two CXO-10 ex- 
posures, the power-law fluxes are (1.2 ± 0.2) x 10 -13 and 
(0.7 ± 0.2) x 10~ 13 erg s" 1 cm" 2 . 

Another possible cause for concern is that tying the 
power-law index between the different observations (except 
for XMM-3) may skew the derived values for the temper- 
ature and non-thermal flux. To gauge to what extent our 
results may be affected by this, we repeated our spectral 
fit twice, with the index fixed at values of 1.3 (similar to 
the value for XMM-3) and 2.5. Fixing the index at 1.3 re- 
sulted in all the temperature values increasing. All but 
two had shifts in the range 0.9-3.8 eV; for CXO-3 and 
CXO-5 the shifts were 8.1 eV and 5.6 eV, respectively. The 
(non-XMM-3) power-law fluxes all decreased by 24%-34%. 
Fixing the index at 2.5 resulted in temperature shifts in 
the range —5.8 to —0.1 eV, except for CXO-3 and CXO-5, 
whose shifts were —16.2 eV and —9.1 eV. The power-law 
fluxes increased by 59%-78%. The derived parameters for 
the cooling curve fits in both these cases are shown in Ta- 
ble |31 no drastic changes are seen. Overall, for most of the 
observations it is unlikely that the individual temperatures 
are skewed by more than ~6 eV due to the tying of the 
power-law index. For CXO-3 and CXO-5, the two observa- 
tions with the largest non-thermal components (excluding 
XMM-3), the effect may be larger. However, we note that 
when the indices for these two observations are allowed to 
vary independently from the others, the temperatures are 
only shifted by -5 eV (CXO-3) and +1.6 eV (CXO-5). 
The values of the two indices in this case are 2.2 ± 0.4 and 
1.8 ± 0.4. The power-law fluxes are likely skewed by less 
than a factor of 2 in all cases. While these possible ef- 
fects on the temperatures and fluxes are unfortunate, the 
situation would probably not be improved by allowing the 
power-law indices to vary independently. In that case the 
indices take on widely varying (and in some cases unphys- 
ically high or low) values, many of which are very poorly 
constrained; almost all the ler confidence intervals for the 
free index values overlap. In addition, the overall fit be- 
comes very unstable. Given all this, it is unlikely that the 
temperatures and fluxes derived with the free indices are 
in general any more accurate than those derived using a 
single "average" value for the index, and for those obser- 
vations where the indices take on extreme values they are 



probably less accurate (although this could be mitigated to 
some extent by constraining the index values to lie in a cer- 
tain physically plausible range, e.g., 1-2.5). We conclude 
that the best solution, although clearly not perfect, is to tie 
the index. 

It is natural to ask whether this possible skewing of the 
temperatures, and/or limitations in our ability to sepa- 
rate the contributions of the two spectral components to 
the total flux, can explain the anomalously high tempera- 
tures seen for XMM-3 and CXO-4. To test this we fitted 
the spectra from these observations individually, fixing the 
temperatures at the value of the best-fit exponential decay 
curve to the main spectral fit (see Figure 2]), and allowing 
the power-law indices to vary. The absorption column was 
fixed at the best-fit value from the main spectral fit. The 
CXO-4 spectrum is well fitted in this manner (xt =0.80 for 
15 dof), although the value of the power-law index is high, 
2.96±0.25 (and very different from the value attained when 
the temperature is also allowed to vary, 0.7 ± 1.1, in which 
case the temperature increases from the original main spec- 
tral fit value). Constraining the temperature at the value 
of the best-fit broken power-law curve (see Figure [5]) re- 
sults in a lower value for the index, 2.55 ± 0.33 (and still 
giving an excellent fit). In light of this, and the fact that 
CXO-4 does not show an anomalously high luminosity or 
non-thermal flux compared to the other observations (in 
contrast to XMM-3), we conclude that it is quite possible 
that the high temperature derived from the main spectral 
fit is in this case simply an artifact of the fitting process, 
and not due to a real increase in the temperature compared 
to CXO-3. Fixing the temperature in XMM-3 to the values 
of the exponential decay or broken power-law curves gives 
unacceptable fits; the xt values are 1.51 and 1.46 for 225 dof 
(P x of 1 x 10~ 6 and 9 x 10~ 6 ). In comparison, the xt value is 
1.18 (P x = 0.033) for the original XMM-3 fit with the tem- 
perature free; much better, but admittedly only marginally 
acceptable. This is partly due to the fact that the model un- 
derestimates the lowest-energy part of the spectrum; this 
can be seen in the pn spectrum in Figure |3J Indeed, al- 
lowing the absorption column to vary leads to a consider- 
ably better fit with xt = 1.11 (P x = 0.13). In this case 
iV H = (1.72 ± 0.05) x 10 22 cm -2 ; the derived temperature 
is 150.2 ±3.2 eV. Untying the absorption column for XMM- 
3 from that of the other observations leads to a significant 
improvement in the overall main fit, from a xt value of 1.12 
for 517 dof (P x = 0.035) to 1.07 (P x = 0.12); the tied iV H 
then assumes a value of 1.99±0.03 cm -2 . Furthermore, do- 
ing this in conjunction with untying the temperatures and 
power-law normalizations within CXO-7 and CXO-10, as 
discussed earlier in this section, leads to a very good fit with 
xl = 1.02 for 508 dof [P x = 0.34). We also note that the 
individual XMM-3 fit, with TVh free and the temperature 
fixed to the value of the best-fit broken power-law curve, 
gives xt — 1-17 (P x = 0.042); a marginally acceptable fit. 
However, the much lower value of the absorption column 
for XMM-3 (in this case 1.59 ± 0.04 cm -2 ), compared to 
the tied value for the rest of the observations, is rather im- 
plausible. Overall, we conclude that the XMM-3 spectra 
are hardly compatible with a temperature in line with the 
overall trend of the rest of the observations, and that there 
is in this case likely a real and significant increase in ther- 
mal flux compared to the previous observation. The high 
luminosity and non-thermal flux also unequivocally show 
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that this observation is quite distinct in behavior from the 
other ones. 

3. DISCUSSION 

We have presented RXTE and Swift observations track- 
ing the final 3 weeks of the 2006-2007 outburst of XTE 
J1701-462, and subsequent Chandra and XMM-Newton ob- 
servations monitoring the source during the first c±800 days 
of quiescence. The transition from active accretion to qui- 
escence is resolved with much better precision than for any 
other cooling NS transient observed after an extended out- 
burst; the end of the outburst is tightly constrained to a 
~4 day window in 2007 August. We fit the spectra ob- 
tained during quiescence with a two-component model con- 
sisting of a NS atmosphere model (thermal component) and 
a power-law model (non-thermal component). The effective 
surface temperature of the NS, derived from the thermal 
component, was seen to decay rapidly in the first ^200 
days of quiescence, which we interpret as the cooling of 
the NS crust toward thermal equilibrium with the core, af- 
ter having been heated by accretion during the outburst. 
Our data set yields a much better sampled cooling curve 
than those of other cooling NS transients observed to date. 
The interpretation of the data, is complicated by an ap- 
parent temporary increase in the temperature ~220 days 
into quiescence. The existence of the non-thermal compo- 
nent in the quiescent spectra also adversely affects how well 
we can constrain the behavior of the thermal component. 
The non-thermal flux from the source has varied irregularly 
throughout the quiescent phase by a factor of ~30, repre- 
senting ~5%-50% of the total flux. Fitting the inferred 
temperatures (excluding the two observations showing the 
temporary increase) with an exponential decay plus a con- 
stant offset yields an e-folding time of ~120^ 3 q days. This 
value is not affected by uncertainties in the NS distance, 
radius, or mass, and uncertainty in the modeling of the 
non-thermal component likely has only a small effect on 
the derived timescale. The exponential fit implies that the 
temperature has reached a roughly constant value of ~125 
eV (assuming the best-estimate distance of 8.8 kpc). Al- 
lowing for uncertainty in the NS distance, mass, and radius, 
and in the modeling of the non-thermal component, gives a 
possible error of ~15 eV. This best-fit baseline temperature 
corresponds to a bolomctric thermal luminosity (redshifted 
and unabsorbed) of ~5.4 x 10 33 erg s _1 . The temperature 
data can also be adequately fitted with a broken power law 
(but not a simple power law) . Taking into account that the 
cooling curve may have flattened out from the power-law 
behavior, the break time is constrained to ~20-150 days 
after the outburst end. 

3.1. Transition to Quiescence 

J 1701 was observed to transition sharply from active ac- 
cretion to quiescence in the final 2 weeks of the outburst. 
The decay rate increased dramatically around 2007 July 
29 and the luminosity subsequently dropped by a factor 
of ^2000 in the final ~13 days before entering quiescence 
(see Figure [5]). The drop is well represented by a sim- 
ple exponential decay with an e-folding time of ~1.7 days. 
This evolution from outburst to quiescence is similar to 
that seen in some other NS X-ray transients, e.g., Aql X-l. 
iCampana et al.l (|1998bl ) analyze observations of Aql X-l at 
the end of a 1997 outburst in which the decay steepened 



suddenly and the luminosity then decreased by 3 orders of 
magnitude in less than 10 days; the decrease is well de- 
scribed by an exponential decay with an e-folding time of 
~1.2 days. This behavior , closely resembl i ng tha t seen for 
J 1701, was interpreted by ICampana et all (|1998bl ) as being 
caused by the onset of the propeller mechanism, which im- 
pedes accretion, and signaling the turning on of a rotation- 
powered pulsar. They interpret the subsequent quiescent 
emission as arising from a shock between the pulsar wind 
and outflowing material from the companion star. How- 
ever, several serious problems with this general interpre- 
tation of the steepening of the outburst dec ay rate in NS 
transie nts have been pointed out (see, e.g., (jonker et al.l 
l2004bl) . These are both observational in nature, such as 
the fact that steepening of the decay has also b een seen in 
black hole transi ents (e.g.. | Jonker et al.ll2004al ). as well as 
theoretical (e.g.. IRappaport et al.ll2004j ). 

3.2. Behavior of the Thermal Component 

It is useful to compare the behavior of J 1701 to that of 
the cooling quasi-persistent transients KS 1731-260 and 
MXB 1659-29. Following their ~12.5 and ~2.5 yr out- 
bursts, the effective surface temperatures of KS 1731-260 
and MXB 1659-29 were seen to decay exponentially with 
e-folding times of 305 ± 47 and 465 ± 25 days, respec- 
tively, reach ing an approximately constan t level in ~1000- 
1500 days (IWiinands et al.l [200ll [20021: IWiinandsl [ 2002; 
I Rutledee et alj|2002t IWiinands et al.l l200a 120041; IWiinandsl 
2004; Cackct t et al.ll2006l 120081) . However, as mentioned in 
Section [TJ a recent observation of KS 1731-260 indicates 
that the sou rce may still be cooling slow ly. It was pointed 
out early on (jWiinands et al.ll2002ll2004D that the KS 1731- 
260 and MXB 1659-29 data indicate d NS crusts with high 
thermal conductivity. More recently, IShternin et al.l (|2007l ) 
compared simulations of deep crustal heating and subse- 
quent crustal cooling to the observations of KS 1731-260, 
and reached the conclusion that low thermal conductivity, 
corresponding to an amorphous (i.e., non-crystalline) crust, 
is inconsistent with the data. This is in agr eement with the 
resul t s of molecular dynamics simula tions (|Horowitz et al.l 
120071 120091; Horowitz fc Berrvl 120091) . which indicate that 
the crust of an accreting NS will form an ordered crys- 
tal. iBrown fc Cummiirg] (|2009t hereafter BC09) construct 
models of the thermal relaxation of a NS crust following 
an extended a c cretio n episode and confirm the finding of 
IShternin et al.l (|2007T > that the thermal conductivity of the 
crust is high. Fitting their models to the observations of 
KS 1731-260 and MXB 1659-29, BC09 are able to place 
constraints on the crust parameters of the NSs; in particu- 
lar, they are able to tightly constrain the so-called impurity 
parameter (which measures the distribution of the nuclide 
charge numbers, Z, in the crust material) for MXB 1659- 
29. They find a low value for the parameter, indicating a 
small amount of impurities in the crust (i.e., material with 
a small spread in Z) and high conductivity. In the case of 
J 1701, the fact that the source seems to have cooled con- 
siderably faster than both KS 1731-260 and MXB 1659-29 
(the derived exponential decay timescale is shorter by fac- 
tors of ~2— 5) strongly indicates a highly conductive crust, 
and possibly suggests low-impurity material. Extracting 
the exact implications of our data requires fitting of the- 
oretical models; however, that is beyond the scope of this 
paper. 

The models of BC09 indicate that the temperature cool- 
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ing curve should actually not follow an exponential decay, 
but should rather be close to a broken power law flatten- 
ing out to a constant at late times. The break is predicted 
to take place a few hundred days after the end of the out- 
burst, and is set by the thermal diffusion time to the surface 
from the depth in the crust at which the material transi- 
tions from a classical to a quantum crystal, close to neu- 
tron drip (which is where neutrons start to leak out of nu- 
clei, thereby producing a neutron gas between the nuclei). 
The break is mainly due to the suppression of the specific 
heat of the crust material with increasing density. Estab- 
lishing whether and when a break in an observed cooling 
curve takes place would provide valuable input for theoret- 
ical models of the NS interior. There is no indication of a 
break having occurred in the cooling curve of KS 1731-260; 
the temperature data are act ually well fitted wit h a single 
power law of slope 0.12±0.01 ([Cackett et al.ll2008| ). A break 
may have occurred for MXB 1659-29 in the first ^500 days 
after the end of the outburst (the fits of BC09 imply a 
break at —300-400 days post-outburst), but the scarcity of 
observations precludes drawing firm conclusions. As can be 
seen in Figure [5j there is a strong indication for a break in 
the cooling curve of J1701 (although see caveats mentioned 
in Section [221) . The time of the possible break, however, 
is only —20-150 days post-outburst, much earlier than the 
break predicted by BC09. It is therefore not clear whether 
the break seen for J1701, if real, is the one predicted by 
BC09, or is perhaps an unrelated extra structure in the 
cooling curve. Judging from Figure 10 in BC09, the time 
of the break can perhaps be decreased to as little as —150 
days by assuming a rock-bottom value of for the impu- 
rity parameter (no impurities in the crust), corresponding 
to a very high thermal conductivity. In addition, BC09 
point out that the timescale of the cooling is proportional 
to R* S M- 2 {1 + z)- 1 . Their model assumes R ns = 11.2 km 
and M ns — 1.62 M Q ; a smaller radius and/or a larger mass 
(although a large mass is somewhat unlikely for the NS 
in J1701; see discussion below) would therefore push the 
break to an earlier time. This suggests that the observed 
break could conceivably be the one predicted by BC09 if 
it is in the upper part of the —20-150 day range. A break 
time in the lower part of the range is not consistent with 
the predicted break. A possible alternative explanation for 
the break is the existence of a strong nuclear heat source 
in the outer crust. Such a source would alter the shape of 
the crust temperature profile in surrounding layers during 
outburst. This could cause a break in the cooling curve 
around a time corresponding to the thermal diffusion time 
from the depth of the source to the surface. Interestingly, 
iHorowitz et al.l ((2)08) calculate that 24 O should fuse at 
densities near 10 11 g cm" 3 , releasing 0.52 MeV per accreted 
nucleon. According to the models of BC09, the thermal dif- 
fusion time to the surface from a depth corresponding to a 
density of 10 11 g cm" 3 is ^100 days (see Figures 5 and 6 in 
their paper); this is consistent with the time of the break 
in the J1701 cooling curve. 

In the models of BC09, the initial slope of the broken 
power-law curve gives a direct measure of the inward flux 
near the top of the crust during outburst. They note that 
observations in the first 2 weeks after the end of an out- 
burst are critical for constraining the depth and strength 
of heat sources in the outermost layers of the crust. J 1701 
was observed 3 times in the first — 16 days of quiescence; 



these observations should be able to provide valuable in- 
put for theoretical models. We note that our data set is 
unique in this respect. MXB 1659-29 and KS 1731-260 
were not observed until at least —3 1 and —48 days, re - 
spectively, had passed in quiescence (jCackett et al.l [20061 ): 
the situation is less clear with EXO 0748-676, since the 
end time of the outburst is on ly constrained to a —7 week 
window (Deg enaar et al.l [2009f ). Our data, in conjunction 
with the results of BC09, allow us to calculate an esti- 
mate for the energy release per accreted nucleon in the 
outermost layers of the crust. Using Equation 12 in BC09 
and our estimate of —0.027 for the pre-break slope of the 
cooling curve gives an outer crust flux during outburst of 
—9.0 x 10 20 erg s _1 cm -2 and, integrating over a surface of 
radius 10 km, a total energy flow rate of ~l.lx 10 34 erg s _1 . 
The observed bolometric energy output of the outburst, 
—1.0 x 10 46 erg (see Section [TTT]) . gives an average lumi- 
nosity of —2.0 x 10 38 erg s _1 over the approximately 19- 
month- long outburst. Assuming an accretion-powered lu- 
minosity, L — eMc 2 with e = 0.2, we get an estimate of 
(M) ~ 1.1 x 10 18 g s" 1 ~ 1.7 x 10" 8 M Q yr" 1 for the av- 
erage mass accretion rate during the outburst. The total 
energy flow rate derived from the early-time slope then cor- 
responds to —11 keV per accreted nucleon. This value has 
a large uncertainty and is probably an underestimate, since 
the accretion rate was lower than the average value in the 
later parts of the outburst. Accounting for the inferred 
outer crust flux therefore requires more energy per accreted 
nucleon in the outermost layers of the crust; these are the 
layers being probed by the cooling right after the end of 
the outburst. Since our data are consistent with a range 
of values for the early-time power-law slope, we note that 
slopes in the range 0.01-0.08 give energies of —4-32 keV 
per accreted nucleon. Heat deposits per nucleon in this 
range are consistent with the energies availa ble from elec- 
tron captures in the outer laye rs of the crust (|Gupta et al.1 
120071: lHaensel fc Zdunik|[2008[ ). We note that most of the 
total heat deposit per nucleon (which is of the order of an 
MeV; see below) is believed to be released in pycnonuclear 
fusion reactions much deeper in the crust, as mentioned in 
Section [1] 

The initial NS surface temperature for J1701 in quies- 
cence seems to be considerably higher than those of the 
other studied cooling sources. For J1701 this temperature is 
-165 eV, compared to -110 eV for KS 1731-260, -130 eV 
for MXB 1659-29, and -120-130 eV for EXO 0748-676, 
as determined by extrapolating the best-fit exponential de- 
cay cooling curve back to the end of the outburst in each 
case. The high temperature for J1701 possibly reflects the 
fact that it accreted at an extraordinarily high level during 
its outburst, much higher than any of the other sources, 
although it could also be related to a higher equilibrium 
temperature for J1701. We note, however, that these values 
could be off by as much as —10-20 eV due to uncertainties 
in the NS distance and radius, and exponential extrapola- 
tion. Moreover, it may not be appropriate to assume expo- 
nential decay at the start of quiescence. If the behavior was 
closer to a power law, the initial temperatures could be sig- 
nificantly higher for KS 1731-260, MXB 1659-29, and EXO 
0748-676, but this is highly uncertain, due to the scarcity 
of early observations. 

It is not clear whether the temperature of J1701 has 
reached its equilibrium value or thermal relaxation between 
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the NS crust and core is still ongoing. The exponential de- 
cay fit to the cooling curve strongly indicates that the tem- 
perature has reached an approximately constant value. The 
situation is less clear when considering the broken power- 
law model (see Figure [5]) . The source could still be in a 
slow decay, or the power-law curve may already have flat- 
tened out as it is expected to eventually do. However, 
such flattening at perhaps —200 days post-outburst would 
have occurred much earlier than for KS 1731-260 and MXB 
1659-29, where fits would indicate a value closer to —1000 
days (BC09), and it is not even clear whether KS 1731-260 
has stopped cooling. Further observations of J1701, 1000- 
2000 days post-outburst, are needed to confirm whether the 
source has reached equilibrium and, if not, to constrain the 
rest of the decay. If the crust has relaxed, then the equilib- 
rium surface temperature of J 1701 is much higher than the 
values reported for KS 1731-260 and MXB 1659-29. The 
best-fit equilibrium temperature for J1701 is 125.0±0.9 eV, 
compared to 70.2 ± 1.2 eV for KS 1731-260 (or perhaps less, 
given the recent observation mentioned in Sectio n [1} and 
54 ± 2 eV for MXB 1659-29 ([Cackett et all [20081 ) . Uncer- 
tainties in the NS distance and radius (as well as, in the 
case of J 1701, uncertainties due to possible effects from the 
non-thermal component on the derived temperatures) make 
these values uncertain by —5-15 eV. This implies a total 
temperature drop for J1701 of —40 eV, similar to that seen 
for KS 1731-260 (not taking into account the most recent 
observation), but considerably less than the —75 eV drop 
for MXB 1659-29 (note, however, our caveat above about 
uncertainties in the initial temperatures). The equilibrium 
temperature is set by the temperature of the NS core, which 
in turn depends on the long-term time-averaged mass ac- 
cretion rate of the NS and the extent to which the core is 
able to cool via neutrino emission. Unfortunately, all the 
information we have on the accretion history of J 1701 is 
the recent outburst and the fact that before the outburst 
the source had probably been in quiescence at least since 
the start of the RXTE mission and its all-sky monitoring 
in 1996 January, although with the ASM we cannot rule 
out long-term activity at luminosities — 10 34 -10 36 erg s" 1 
or short-term activity at higher luminosities. The recur- 
rence time for outbursts in this system is therefore un- 
known, and we do not know whether the observed outburst 
is representative of typical behavior for the source. The 
possible equilibrium temperature of J1701 corresponds to a 
bolometric thermal luminosity of —5.4 x 10 33 erg s" 1 . This 
would be among the highest luminosities seen for a quies- 
cent NS-LMXB, si milar to those seen for A ql X-l and 4U 
1608-52 (see, e.g.. iHeinke et al.ll2007l 120091 and references 
therein) . 

To get some indication of how our results for J 1701 com- 
pare with theoretical predictions of the quiescent thermal 
luminosities of accreti ng NS transients, we use the results of 
lYakovlev et al.1 ([20041 ) . They compute the quiescent bolo- 
metric thermal luminosity as a function of long-term time- 
averaged mass accretion rate for several models of accreting 
NSs warmed by deep crustal heating. They do this for sev- 
eral different equations of state, for NS masses between 1.1 
Mq and —2.0 M , and for two models of a heavy-element 
accreted envelope: nuclear burning ashes composed of 56 Fe 
with a total heat deposit of 1.45 MeV per accreted nucleon 
(jHaensel fc ZdunikUl990T ). and 1 TO Pd ashes with a heat de - 
posit of 1.12 MeV per nucleon (jHaensel fc Zdunikl [20031) . 
They also take into account a possible He layer on top of the 



heavy - element envelope. We note that Ua ensel fc Zd unik 
(2008) have since revised the estimates of the heat deposits 
to -1.9 MeV for 56 Fe and -1.5 MeV for 106 Pd. We com- 
pare our results to the computed luminosity curves for a 
NS of mass 1.1 Mq (the authors note that the curves are 
nearly identical for 1.3 Mq, and are insensitive to the as- 
sumed equation of state among those they consider), and 
for a heat deposit of 1.45 MeV per nucleon (see Figure 5 in 
lYakovlev et al.ll2~004l ). Given the quoted 15% uncertainty 
in the best-estimate distance to J1701 and various uncer- 
tainties arising from, e.g., the assumed values of the NS pa- 
rameters, the choice and fitting of the spectral model, and 
the extrapolation of the thermal component outside the ob- 
served energy range (see Sections [2.51 and [2T6]) . we consider 
a representative range of (3-9) x 10 33 erg s _1 for our as- 
sumed equilibrium bolometric thermal luminosity. For the 
upper limit luminosity of 9 x 10 33 erg s _1 and no He layer, 
we get an upper limit accretion rate of— 2xl0 -9 Mq yr . 
For the lower limit luminosity of 3 x 10 33 erg s _1 and the 
case of a thick He layer, we get a lower limit accretion rate 
of —4 x 10 -11 Mq yr -1 . This range is consistent with con- 
strai nts on the accretion rates of many NS transients (see, 
e.g., IHeinke et all 120071 120091 ) . Assuming that the 2006- 
2007 outburst is typical for J 1701, these long-term averaged 
accretion rates, in conjunction with our estimated average 
outburst rate (see discussion earlier in this section), give 
recurrence times for outbursts in the system ranging from 
— 10 years to —700 years. Our data for J1701 are therefore 
entirely consistent with standard cooling in a low- mass NS, 
and there is no need to assume a higher-mass NS with en- 
hanced cooling. In fact, significantly enhanced cooling is 
somewhat unlikely, since in that case a very high long-term 
average accretion rate would be required to keep the NS as 
warm as it is observed to be; this suggests a rather low- 
mass NS in J 1701, especially for NS models assuming nu- 
cleon or nucleon-hyp eron matter in the core (see Figure 3 in 
Yako vlev et~all 120041 and Figure 12 in lYakovlev fc Pethickl 
20041) . We note that our estimate of the outburst accretion 
rate could well be off by a factor of —5-10 given uncertain- 
ties in the energy output of the outburst and the radiative 
efficiency of the accretion, and substantial uncertainties are 
of course associated with the theoretical calculations. Nev- 
ertheless, we conclude that if the crust of J1701 has indeed 
already reached (or is close to reaching) thermal equilib- 
rium with the core, then a rather low-mass NS with a near- 
standard cooling scenario is more likely than a high-mass 
star with significantly enhanced cooling. 

3.3. Behavior of the Non-thermal Component 

We speculate that the increase in temperature seen for 
the third XMM-Newton observation (XMM-3) is due to an 
additional spurt of accretion. This is supported by the 
fact that XMM-3 has a large non-thermal flux, which is 
suggestive of ongoing accretion given the fact that promi- 
nent non-thermal components are commonly seen in spec- 
tra from accreting NS systems at low luminosity (see, e.g., 
iDi Salvo fc STeliar[200l iBarret et alJl2000h . During such 
an accretion spurt the surface of the NS may have been 
subjected to some shallow and temporary reheating. Fur- 
thermore, [ZamjjierreFal. (1995) show that low-level accre- 
tion onto a NS surface can produce radiation with a ther- 
mal (hardened blackbody-like) spectrum; additional ther- 
mal flux of this sort could also explain the increase in the 
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inferred effective temperature. The subsequent Chandra 
observation (CXO-4) seems to have an anomalously high 
temperature as well, but has a much smaller non-thermal 
flux than XMM-3. As discussed in Section [2. 6.11 the high 
effective temperature derived for CXO-4 may be a result of 
the non-thermal component in the spectrum skewing our 
estimate of the thermal flux; a temperature in line with 
the overall observed decay is also consistent with the data. 
However, if the heightened thermal flux in CXO-4 is real 
and due to accretion, then the fact that this observation 
does not show a large thermal flux compared to the other 
observations (in contrast to XMM-3) is possibly indicative 
of some change in the properties of the accretion flow be- 
tween XMM-3 and CXO-4. In any case, the substantial 
decrease in both the thermal and non-thermal flux com- 
pared to XMM-3 would point to a significant decrease in 
the rate of this possible accretion. We note that no activity 
was seen with the RXTE ASM around the time of XMM-3 
(or at any other time since the start of quiescence). This 
is hardly surprising, since an increase in luminosity of sev- 
eral orders of magnitude would be required for an ASM 
detection. 

In general, a high degree of variability is seen in the non- 
thermal flux from J 1701. Variability by a factor of ~30 
(~10 when excluding XMM-3) is observed between the in- 
dividual observations obtained throughout the first ~800 
days of the quiescent phase. Indications of possible vari- 
ability on a timescale of ~1— 3 days is seen in the data 
from CXO-7 and CXO-10. As irregular variability on vari- 
ous timescales is a common characteristic of accretion, this 
points to accretion of some sort as a strong candidate for 
the cause of the non-thermal component. We searched for 
short-term variability (within observations) by constructing 
light curves for the XMM-Newton observations, but did not 
see indications of variability. For each Chandra observation 
we performed a one-sided K-S test comparing photon ar- 
rival times to a uniform count rate model. This did not 
reveal any evidence for variability at a confidence above 
90%, with the exception of the first of the two exposures 
in CXO-10, which showed variability at a 97.7% confidence 
level (discussed in Section [2 .6. 1[) . However, due to the low 
count rates in our observations this does not set strong con- 
straints on the possible existence of short-term variability. 

There seems to be a significant difference between the 
value of the power-law index seen for XMM-3 (~1. 3-1.4) 
and those seen for the Swift and RXTE observations dur- 
ing the transition to quiescence (~1.7 and ~1.9 for the first 
and second Swift observations, and even higher for the third 
one, although highly uncertain due to very few counts; the 
RXTE observations almost all had index values in the range 
~1.8-2.5). Given that the outburst was still ongoing when 
the Swift and RXTE observations were made, it is reason- 
able to assume that those spectra are dominated by flux 
due to accretion. This may indicate a different origin for 
the non-thermal component, or a difference in the accretion 
flow compared to the outburst, for XMM-3. We do note, 
however, that the luminosity in XMM-3 was lower, by a 
factor of 2 or more, than the luminosity in any of the Swift 
and RXTE observations; the validity of a direct compari- 
son is therefore questionable. The value of the combined 
power-law index for the other observations in quiescence is 
in line with those from the Swift and RXTE observations, 
but this combined index does not give us information about 
the index values for individual observations. The fact that 



the combined index value seems to be different from the 
XMM-3 one may also point to some difference in nature 
between the non-thermal component in XMM-3 and that 
in the other quiescent obse rvations. 

A study by Jonkcr ct al. (2004b) suggested that the frac- 
tional contribution of the non-thermal flux to the total un- 
absorbed 0.5-10 keV flux in non-pulsing NS-LMXB tran- 
sients evolves as a function of total 0.5-10 keV luminos- 
ity. Ba sed on data from several NS-LMXBs, Uonker et al.l 
(2004b) suggest that the fraction decreases with decreas- 
ing luminosity down to a minimum of 10%-20% around 
r^2 x 10 33 erg s _1 , and increases again at lower luminosi- 
ties to values of 60%-70% (see Figure 5 in their paper). In 
Figure [6j we show the fractional contribution of the power- 
law component to the total 0.5-10 keV flux as a function 
of the total luminosity in the same band. The power-law 
fraction varies rather randomly between ~5% and ~50%, 
an d does not seem to ex hibit behavior of the sort seen in 
the Uonker et al.l (|2004b[) study. 

4. SUMMARY 

We have presented RXTE and Swift observations tracking 
the final 3 weeks of the 2006-2007 outburst of the super- 
Eddington NS transient XTE J 170 1-462, as well as Chan- 
dra and XMM-Newton observations covering the first ~800 
days of the subsequent quiescent phase. The source tran- 
sitioned sharply from active accretion to quiescence, with 
the luminosity decreasing by a factor of ^2000 in ~13 days. 
The end of the outburst is tightly constrained to a ~4 day 
window in 2007 August. 

We fitted the Chandra and XMM-Newton spectra with a 
two-component model consisting of a NS atmosphere model 
(thermal component, interpreted as being due to emission 
from the NS surface) and a power-law model (non-thermal 
component whose origin is uncertain). The effective surface 
temperature of the NS, inferred from the thermal compo- 
nent, was seen to decay rapidly in the first ~200 days of 
quiescence. We interpret this as the NS crust cooling after 
having been heated and brought out of thermal equilibrium 
with the core during the outburst. The interpretation of 
the data is complicated by an increase in the derived tem- 
perature ~220 days into quiescence. The existence of the 
non-thermal component also adversely affects our ability to 
constrain the thermal component. 

Fitting the derived temperatures with an exponential de- 
cay cooling curve plus a constant offset (excluding the two 
observations affected by the apparent temperature increase) 
we derive an e-folding time of ~ 120^2° days with a best- fit 
offset of ~125 eV. This baseline temperature is uncertain 
to ~15 eV due to uncertainties in the distance, radius, and 
mass of the NS, and in the modeling of the non-thermal 
component. The short decay timescale strongly indicates 
high thermal conductivity in t he NS crust, and possibl y 
suggests low-impurity material (Brown & dimming 2009) . 

The temperature data can also be well fitted with 
a more physically motiv ated broken power-law model 
(|Brown fc Cumminel [2009) . The data show a strong in- 
dication for a break in the power law; the time of the break 
is ~20-150 days after the end of the outburst. This is 
considerably earlier than the break predicted by theory, 
which is mainly due to a change in heat capacity in the NS 
crust where the material transitions from a classical to a 
quantum crystal (|Brown fc Cum ming 20Q9|). The observed 
break may therefore have a different origin; a possible al- 
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ternative explanation is a strong nu clear heating source i n 
the crust, e.g., the fusion of 24 (|Horowitz et alj [2008h ■ 
The initial slope of the power-law cooling curve is a direct 
measu re of the inward flux in th e outer crust during out- 
burst (jBrown fc Cummind I2009T) . Our measured slope, in 
conjunction with an estimate for the accretion rate dur- 
ing the outburst, yields an estimate for the heat deposit 
per accreted nucleon in the outermost layers of the crust 
which is consistent with theoretical pre dictions for the en- 
ergy available from elec tron captures (|Gupta et al.l [2007; 
lHaensel fc Zdunikl 120081 ). 

Further observations are needed to determine whether the 
crust is still cooling slowly or has already reached thermal 
equilibrium with the core at a surface temperature of ~125 
eV. The latter would imply an equilibrium bolometric ther- 
mal luminosity of ~5 x 10 33 erg s" 1 for an assumed distance 
of 8.8 kpc. This would be among the highest quiescent ther- 
mal luminosities seen from a NS-LMXB, and may indicate a 
rather low-mass NS without significantly enhanced cooling. 

The non-thermal component has varied irregularly 
throughout the quiescent phase by a factor of ~30; indi- 
cations of possible variability have been seen on timescales 
as short as ~1 day. The fractional contribution of this 
component to the total flux has been between ~5% and 
~50%. We speculate that the non-thermal component in 
XTE J1701-462 arises from residual accretion, and that the 
increase in derived temperature ~220 days into quiescence 
(which was accompanied by a large increase in non-thermal 
flux) was due to a spurt of increased accretion, possibly 
causing some shallow and temporary reheating of the NS 
surface and/or releasi ng radiation with a thermal spectrum 
(jZampieri et alj|1995j ). 

The observed behavior of XTE J1701-462 during the 
post-outburst quiescence seems to be quite different from 
that of other cooling NS transients observed after extended 
outbursts: the derived exponential decay timescale of the 
effective surface temperature is much shorter; both the ini- 
tial and final surface temperatures are likely significantly 
higher; there is a strong indication for a break in the cool- 
ing curve when fitting the temperature data with a (bro- 
ken) power law; a significant temporary increase in thermal 
flux was observed after more than 200 days of quiescence; 
finally, the source has exhibited a prominent non-thermal 
component in its spectrum throughout the quiescent phase. 
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